PSR J2021+3651: A NEW 7-RAY PULSAR CANDIDATE 
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ABSTRACT 

The COS B high, energy 7-ray source 2CG 075+00, also known as GeV J2020+3658 or 3EG J2021+3716, has avoided 
identification with a low energy counterpart for over twenty years. We present a likely identification with the discovery 
and subsequent timing of a young and energetic pulsar, PSR J202 1+3651, with the Wideband Arecibo Pulsar Processor 
at the Arecibo Observatory. PSR J202 1+3651 has a rotation period P ~ 104 ms and P ~ 9.6 x 10^^^, implying a 
characteristic age Tc ~ 17 kyr and a spin-down luminosity E ~ 3.4 x 10^^ erg s~^. The pulsar is also coincident 
with the ASCA source AX J2021. 1+3651. The implied luminosity of the associated X-ray source suggests the X- 
ray emission is dominated by a pulsar wind nebula unresolved by ASCA. The pulsar's unexpectedly high dispersion 
measure (DM ~ 371 pc cm~^) and the d > 10 kpc DM distance pose a new question: is PSR J2021+3651 an 
extremely efficient 7-ray pulsar at the edge of the Galaxy? This is a question for AGILE and GLAST to answer. 
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INTRODUCTION 

Finding low energy counterparts to 7-ray sources is difficult because they have large positional uncertainties, 
typically ~ 1° acros s. Consequently, the m ajority of the high energy 7-ray sources observed by COS B and EGRET 
remain unidentified ( Hartman et al. . 19991) . Young pulsars are viable candidates for many of these sources, and they 
remain the only discrete Galactic source class (other than the Sun) unambiguously shown to emit radiation in the 
100 MeV - 10 GeV range. The increased sensitivity to pulsars with high dispersion measure (DM), provided by 
new pulsar backends such as the Parkes multibeam system (Mancheste r et al., 2001), has led to a number of dis 



coveries recently of young pulsars coincident with known 7-ray sources (iD'Amico et al.L l200ll: ICamilo et al.L 12001 



iHaloern et al 



(iCamilo et al 



200?). The recent detection of a young, low luminosity, radio pulsar in the supernova remnant 3C58 



20021) suggests many more faint radio pulsars await discovery in deep, targeted, searches. 



The COS B source 2CG 075+00 has long been considered a pulsar candidate due to its hard spectru m. EGRET re- 



solved it into two sources. The fainter one, 3EG J2016+3657, was identified as a probable blazar bv lMukherjee et. al 
( 20001) . A ROSAT image and an ASCA GIS image based on the second EGRET catalog position of the brighter 
source, 3EG J2021+3716, failed to yield any obvious candidates. These sources are in the crowded Cygnus re- 
gion and the likelihood analysis used to find ^-ra y source positions is sensitive to the assumed number of surround- 
ing sources. 
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Romani. and Kawail ( 2001 ) (hereafter, RRK) noted that the catalog of sources above 1 GeV 



dLamb and MacornbT 1997h contained two sources in this region that were not in the list of 3 o" sources, based on 



100 MeV and above maps, used to construct the 3rd ^G/^fT catalog. They rederived the positional contours of 3EG 
J2021+3716 using 1 GeV and above maps and all of the known nearby sources in their analysis. An ASCA image of 
the new position yielded two hard point-like sources (Figure 0, embedded in softer extended emission. One is iden- 



Table 1 . Measured and derived parameters for PSR J2021 +3651 



Parameter Value 



Right ascension a (J2000) 


20^ 21"^05^214(51) 


Declination 8 (J2000) 


+36° 51'08".42(72) 


Galactic longitude / 


75.°23 


Galactic latitude h 


+0.°11 


Pulse period P (s) 


0.103722225646(12) 


Period derivative P 


9.56118(60) X 10-14 


Pulse frequency v (s^^) 


9.6411351933(11) 


Frequency derivative z> (s^^) 


-8.88726(54) x lO'^^ 


RMS residual 


292/iS 


Epoch (MJD) 


52407.389 


Dispersion Measure DM (pc cm ^) 


369.12(58) 


Pulse width at 50% of peak W50 (ms) 


9.9 


Pulse width at 10% of peak wio (ms) 


18 


Flux density at 1425 MHz (mJy) 


~ 0.1 


Spin-down luminosity E'^ (erg s^^) 


3.4 X 10^6 


Surface dipole magnetic field (G) 


3.2 X 10^2 


Characteristic Age = \P/P (kyr) 


17 


" E = Att^IP/P^ with / = g cm^ 




* Assuming standard magnetic dipole spindown: B 


= 3.2 X 10i9(PP)i/2 


Gauss (Manchester and Tavloi. 1977) 





tified with the Wolf-Rayet + 0-star binary system WR 141. The second, AX J202 1.1+3651, is moderately absorbed 
(riH = 5 X 10^1 cm^^) with a hard, power-law spectrum (photon index T = 1.7). 

By targetting the X-ray source AX J20 21.1+3651 with the A recibo radio telescope, we have discovered the 
young and energetic pulsar PSR J202 1+3651 (^Rob erts et al. . 2002ah . Here we present a timing solution of the pulsar 



based on eight months of monitoring with Arecibo, and discuss the likely prospect that this is the counterpart to 
2CG 075+00, which has eluded identification for over twenty years. 

OBSERVATIONS AND ANALYSIS 

PSR J2021+3651 was observed 17 times between MJD 52305 and 52545 using the Wideband Arecibo Pulsar 
Processor (WAPP) at the Arecibo Observatory. The WAPP is a fast-dump digital correlator with adjustable bandwidth 
(50 or 100 MHz) and variable numbers of lags and sample times (for details see Dowd, Sisk, and Hagen 2000). 
Our observations were made at 1.4 GHz with 100 MHz of bandwidth, 512 lags, 200 fis sampling, and summed 
polarizations. The 16-bit samples were written to a disk array and then transfered to magnetic tape for later analysis. 

Analysis of the observations was done using the PRESTO software suite ( Ransom. .200 l.V . Integrated pulse 
profiles (Figure IJJi from our observations were convolved with a template profile to extract 69 topocentric times of 
arrival (TOA). Using TEMPO^, the topocentric TOAs were converted to TOAs at the solar system barycenter at infinite 
frequency and fit simultaneously for pulsar period, period derivative, DM, RA, and Dec with a residu al rms of 292 us . 
Table^gives the measured and derived parameters for PSR J2021+3651, and is updated from Roberts et al. ( 2002ah . 



The quoted errors are three times those given by TEMPO in order to compensate for systematic errors otherwise 
unacounted for. 



^http://www.physics.mcgill.ca/~ransom/ 
^See http://pulsar.princeton.ed u/tempoi 
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Fig. 1 . ASCA GIS 2-1 keV images of tine COS B region 2CG 075+00 superimposed on an Einstein IPC 
0.1-4.5 keV image. Tine liglit grey contours are tine 68%, 95%, and 99% confidence regions of tine tliird 
SG/e^r catalog sources 3EG J2021 +371 6 (left) and 3EG J201 6+3657 (right). The dark contours are for the 
sourc e GeV J2020+3658 and are derived from the > 1 GeV photons only (Roberts. Romani. and Kawai, 



20011) . The oval is the positional error derived for GeV J2020+3658 by L amb and Macomb. 11997.) . The 
circle centered on AX J2021 .1+3651 indicates the size of the 3' Arecibo beam. 



DISCUSSION 



The DM of PSR J2021+365 1 is by far the highest known in the Galactic long itude range 55° < / < 80°, which 
is mai nly an inter-spiral arm direction. Using the recent Cordes an d Laziol ( 2002 ) update to the Tavlor and Corded 



(119931) dispersion measure model gives a distance of ~ 12.4 kpc, at the outer edge of the last spiral arm used in the 
model. It is possible that there are further contributions from clouds in the Cygnus region not included in the model, 
where there is known to be excess gas at d 1.5 kpc (J. Cordes, private communication, 2002). However, there are 
no obvious HII regions within the Arecibo beam seen in either Very Large Array (VLA) 20-cm radio or Midcourse 
Space Experiment (MSX) 8.3 //m images (available from the NASA/IPAC Infrared Science Archive). 

The likelihood that PSR J202 1+3651 is unrelated to the X-ray source is extremely small given the rarity of 
such young, energetic, pulsars and the consistency of the timing position with the ~ 1' positional uncertainty 
of the ASCA source. However, the high DM is somewhat surprising given the X-ray absorption quoted by RRK 
(nH=(5.0±0.25) X 10^^ cm~^), where the errors represent the 90% confidence regio n. The total Galactic HI colu mn 
density in this direction as estimated from the FTOOL nh, which uses the HI map of iDickey and Lockmanl(ll99nh . 
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Fig. 2. 1.4 GHz pulse profile for PSR J2021+3651 adopted from iRoberts et al.l j2002ah . The error bar 
represents the 1 a uncertainty. 

1.2 X 10^^ cm~^. This should be a good approximation if the source is truly at the far edge of the outer spiral arm. 
Noting that the ASCA image shows faint, softer emission in the region (Figure Q, and given the likely possibility of 
either associated thermal X-ray flux from a supernova remnant or a nearby massive star, we fit the ASCA spectrum 
of RRK, adding a thermal component to t he absorbed power-law model. Accountin g for ~ 4% of the photon flux 
with a MEKAL thermal plasma mo del (see 'Liedahl. Osterh eld. and GoldsteinI ( '^^) and references therein) of tem- 
perature kT ~ 0.1 keV in XSPEC (lAmaud^.l99d) statistically improves the fit (F-test chance probability of 2.5%). 
The best-fit absorption for this three component model is nH=7.6xlO^^ cm~^ with a 90% confidence region of (4.1 
- 12.3) X 10^^ cm^^, consistent with the total Galactic column density. The best-fit photon index is F = 1.86, still 
consistent with the 1.47 - 2.01 range in RRK derived from the simple absorbed power-law model. Hence the X-ray 
absorption does not force us to adopt a smaller distance than is suggested by the DM and the value of nn may be 
contaminated by thermal emission. Upcoming Chandra observations of the source may well clarify this. 

For a distance dio = d/lQ kpc, the inferred isotropic X-ray luminosity Lx = 4.8 x lO'^^dfo erg (2 - 
10 keV). The X-ray efficiency r/x = Ly^/ E is O.Oldfg- Compared to the total pulsar plus nebula X-ray lu r ninosity of 



other spin-powered pulsars, this is somewhat high, but within the observed scatter (iPossenti et al.L 120021 : IChevalier 



ioool) . Upcoming Chandra observations will determine what fraction of the X-ray flux comes from a compact nebula, 



and will allow us to search for X-ray pulsations. 

2CG 075-1-00 = 3EG J202 1-1-37 16 = GeV J2020H-3658 is the 10th brightest Galactic source above 1 GeV. It has a 
hard (photon index F = 1.86 ± 0.10) spectrum with no sign of a break out to 10 GeV, and has low variability, similar 
to the known 7-ray pulsars. PSR J202 1-1-36 5 1's position is outside the 99% contour of the 3EG catalog position, but 
consistent with the 95% contour of both the Lamb and MacombI ( 1997 ) GeV catalog position and the RRK position. 



The only other potential counterpart within any of the error contours is WR 141. Although such binary systems 
have the potential for producing 7-rays from colliding winds, no convincing association with a known high-energy 
(> 100 MeV) 7-ray source has yet been made for a Wolf-Rayet system (Romero. Benaglia. and Torres, 1999j- The 
pulsar's positional coincidence coupled with the high inferred spin-down luminosity strongly suggests this pulsar 
emits 7-rays. Unfortunately, confirming this by folding archival EGRET data is problematic due to the likelihood of 
significant past timing noise and glitches, which make the back-extrapolation of the rotational ephemeris uncertain. 
Assuming the distances given by the .Cordes and Lazia (.2002) model, the known 7-ray pulsars have 7-ray ef- 



ficiencies t]^ = Lj/E mostly between 0.0001 and 0.03 (assuming 1 sr beaming). Assuming the 7-rays are 100% 
pulsed, the inferred 7-ray efficiency for PSR J2021+3651, r/^ = O.lSdfo in t he 100 MeV to 10 GeV range, would 
be by far the most efficient 7-ray pulsar. Note h owever that PSR B 1 055 — 52 ( Thompson et al. , 19991) has a similar 
efficie ncy if one uses the DM distance given by Tavlor and Cordes (1993), instead of the newer ICordes and Laziol 
(l2002h model. PSR B 1706— 44, which has nearly identical spin parameters, has a 7-ray efficiency of ~ 0.01, while 
Vela with a similar spin-down energy has a low efficiency of 0.0004. Even if the pulsar is located within the 
Perseus arm at a distance of 5 kpc, it would still be very efficient. If it is the 7-ray source, it would further strain the 



relationship of efficiency to spin-down energy r/^ oc L_,J^'^ expected from theory ( Zhang and Hardin'3. 2000h . While 



there is currently no observational evidence for a closer distance, increased DM from an intervening source in this 
relatively crowded direction would not be surprising. We note that the DM derived distances for several other young 
pulsars rece ntly di scovered within EGRET error boxes, if the true counterparts, also tend to have high inferred 7-ray 
efficiencies (Roberts I (l2002bh . 

We plan to continue timing observations of this source, which will refine the position, provide a contemporaneous 
timing solution for future high-energy pulse searches, and allow us to monitor for glitches. AGILE and GLAST should 
be able to unambiguously determine if the pulsar is the counterpart of the 7-ray source. 
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